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ABSTRACT 

We present results of five years of optical (UBVRI) observations of the very- 
high-energy gamma-ray blazar 1ES 1011+496 at the MDM Observatory. We cali- 
brated UBVRI magnitudes of five comparison stars in the field of the object. Most 
of our observations were done during moderately faint states of 1ES 1011+496 
with R > 15.0. The light curves exhibit moderate, closely correlated variabil- 
ity in all optical wavebands on time scales of a few days. A cross-correlation 
analysis between optical bands does not show significant evidence for time lags. 
We find a positive correlation (Pearson's r = 0.57; probability of non-correlation 
P(> r) ~ 4 x 10~ 8 ) between the R-band magnitude and the B - R color in- 
dex, indicating a bluer-when-brighter trend. Snap-shot optical spectral energy 
distributions (SEDs) exhibit a peak within the optical regime, typically between 
the V and B bands. We find a strong (r = 0.78; probability of non- correlation 
P{> r) rj 10~ 15 ) positive correlation between the vF v peak flux and the peak 
frequency, best fit by a relation uF^ k oc u^ k with k = 2.05 ± 0.17. Such a cor- 
relation is consistent with the optical (synchrotron) variability of 1ES 1011+496 
being primarily driven by changes in the magnetic field. 

Subject headings: galaxies: active - - BL Lacertae objects: individual (1ES 
1011+496) — radiation mechanisms: non-thermal 
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Introduction 



Blazars are the most violent class of active galactic nuclei, consisting of flat-spectrum 
radio quasars (FSRQs) and BL Lac objects (named after their historical prototyp e, BL Lacer- 
tae). They exhibit rapid v ariability down to time scales as short as a few minutes (lAharonian et al 
20071 ; lAlbert et al.l l2007al ) . Their observed flux is dominated by a non-thermal continuum 
exhibiting two broad spectral bumps: A low-frequency bump from radio to UV - X-ray 
frequencies, and a high-frequency component from X-ray to 7-rays. In the framework of 
relativistic jet models, the low-frequency (radio - optical/UV) emission from blazars is in- 
terpreted as synchrotron emission from nonthermal electrons in a relativistic jet. The high- 
frequency (X-ray - 7-ray) emission could either be produced via Compton upscattering of 
low frequency radiation by the same electrons responsible for the synchrotron emission (lep- 



tonic jet models; for a recent review see, e.g., iBottcherl 120071 ). or due to hadronic processes 



initiated by relativistic p rotons co-accelerated with the electrons (had ronic models, for a 



recent discussion see, e.g.. iMiicke Protherodl2001t IMiicke et al.ll2003l ) 



To date, about 30 blazars have been detected in very high energy (VHE, > 100 GeV) 7- 
rays with ground-based air Cerenkov telescope facilitie^J. Most of these TeV blazars belong 
to the sub-class of high-frequency peaked BL Lac objects (HBLs). They are characterized by 
a synchrotron spectrum peaking at frequencies > 10 15 Hz, i.e., in the UV or X-ray range, 



and 7 -ray peaks at z/P k > 10 25 Hz, i.e., typically beyond the Fermi energy range ( lAbdo et al. 



20101 ) . VHE 7-rays from sources at cosmological distances can b e absorbed by the extra- 



galactic background light (EBL) due to 77 pair production (e.g . . iDwek fc Krennrichl 12005 



Stecker et al. 2006; Franceschini et al. 2008 



Gilmore et al.l 120091 : iFinke et al.ll2010h . Hence 



the EBL absorption of distant VHE 7-ray sources may provide a probe of the spectrum and 
cosmological evolution of the EBL, which is notoriously difficult to measure directly due to 
bright foregrounds. 

The BL La c object 1ES 1011+4 96 was detected as a VHE 7-ray emitter by MAGIC in the 
spring of 2007 ( lAlbert et al.ll2007bl ). Follow-up optical spectroscopy at the MMT confirme d 
the previously uncertain redshift of z = 0.212 ± 0.002 for this source (lAlbert et al.ll2007bl ). 
At the time of its VHE detection, 1ES 1011+496 was the most distant VHE 7-ray source 
known with a well-determined redshift, and to date it still ranks among the top five. It might 
therefore offer a prime opportunity for studying the EBL through its absorption signature at 
VHE 7-rays. However, in order to exploit this opportunity, a thorough understanding of the 
intrinsic spectral energy distribution (SED) of the source, constrained through observations 
at lower (radio through GeV 7-ray) frequencies, is essential. For this purpose, we present 



1 For a complete list of VHE 7-ray sources see http://tevcat.uchicago.edu 
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here a detailed study of the optical spectral variability of this object with data gathered over 
the course of 5 years at the 1.3m McGraw-Hill Telescope of the MDM Observatory on Kitt 
Peak, Arizona. 

The VHE 7-ray detecti on of 1ES 1011+496 was triggered by a large optical outburst 
of the source in 2007 March ([Albert et al.ll2007bl ). This object is regularly monitored in the 
optical R-band through the Turku blazar monitoring program led by Kari Nilsson, with the 
1.03 m telescope of the Tuorla Observatory in Finland, as well as the 35 cm telescope of the 
KVA Observatory on La Palma, Canary Islands, Spairo Apart from this, 1ES 1011+496 has 
so far received rather little attention by optical observers, and its optical spectra l variability 
has r emained unexplored. The object has been observ ed in X-rays by Einstein ( jElvis et al. 
19921 ) and more recently, Swift/XRT (lAbdo et al.ll2010l ). In all observations, it shows a steep 
X-ray spectrum, indicating the dominance of synchrotron emission in the X-ray regime. 
1ES 1011+496 may also be associated with the EGRET 7-rav source 3EG J1009+485 5 
( jHartman et al.lll999l ). although this association is uncertain (ISowards-Emmerd et al.ll2003l ). 
The object is clearly detected by Fermi and listed in the Fermi 3-month catalogue as the 
source 0FGL J1015. 2+4927. The Fermi data rev eal a rising vF u sp ectrum (i.e., photon index 
T < 2) in the 100 MeV - 30 GeV energy range JAbdo et alihoioh . 



Optical/UV observations were performed by Swift/UVOT in May 2008 ( lAbdo et al. 



20101 ). during the rising phase of an optical outburst similar to the one triggering the MAGIC 
discovery observations in 2007. During those observations, Swift/UVOT measured a rising 
vF v optical/UV continuum spectrum. This, together with the hard Fermi spectrum, justified 
the classification of 1ES 1011+496 as an HBL. 

Over the course of ~ 5 years, we have collected optical multi-band (UBVRI) photometric 
data on 1ES 1011+496. As there are no comparison stars with reliably calibrated UBVRI 
magnitudes in the field of view of this object, we calibrated the magnitudes for 5 comparison 
stars. In £j2]we describe the observations, data reduction, the calibration of comparison-star 
magnitudes, as well as general features of the light curves. Results of a spectral-variability 
study are presented in ^J3j We performed a cross-correlation analysis between the variability 
patterns in different optical bands, which we describe in §H We summarize and discuss our 
results in 



2 Daily light curves are posted at |http : //users . utu . f i/kani/ lm/ index . html 
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Fig. 1. — Finding chart (R band) of the field around 1ES 1011+496 with 5 comparison stars 
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2. Observations, data reduction, and light curves 

Optical (UBVRI) data were collected at the 1.3m McGraw-Hill Telescope of the MDM 
Observatory on the south-west ridge of Kitt Peak, Arizona, during 9 ~ 1-week observing 
runs in 2005 April, September, and December, 2008 April and May, 2009 February and 
May, and 2010 April and July. The telescope is equipped with standard Johnson- Cousins 
UBVRI filters. Exposure times for science frames on 1ES 1011+496 were between 60 - 
180 s, depending on filter and atmospheric conditions. All frames were bias-subtracted and 
flat-field corrected using standard routines in IRAF. 

In order to be able to perform relative photometry with comparison stars in the field 
of view, well-calibrated UBVRI magnitudes of those stars need to be known. Since no such 
information is currently available in the literature, we calibrated UBVRI magnitudes of five 
comparison stars in the field of view, as indicated in Fig. [TJ For the calibration, we obtained 
a total of 30 (2 per filter per standard) exposures of the Landoldt Equatorial standards 



flLandoldtlll992D PG 1525-071, PG 1633+099, PG 1657+078 on the night of April 28, 2005, 
which provided photometric observing conditions. Instrumental magnitudes of these three 
Equatorial standards as well as our comparison stars in the field of 1ES 1011+496 were 
extracted using the phot routine within the IR AF package DAOPHOT . Following standard 



procedures for IRAF photometric calibrations ( iMassey fc Davislll992l ). we used the routine 
fitparams to solve the transformation equations to evaluate the calibrated, physical magni- 
tudes of our standard stars. The resulting calibrated magnitudes are listed in Table [TJ 

After calibration of our comparison stars, we extracted instrumental magnitudes of the 
comparison stars and 1ES 1011+496 using the phot routine within the DAOPHOT package 
of IRAF, and converted instrumental to physical magnitudes assuming that the difference 
between instrumental and physical magnitudes is the same for the object and all comparison 
stars. The resulting light curves for all observing runs combined are displayed in Fig. [2j 
The figure illustrates that for most of our runs, the object was in a rather faint optical state 
with R > 15.0, and shows very moderate variability. 

Table 1. Calibrated Magnitudes of Comparison Stars in Fig. [TJ 



Star U 


B 


V 


R 


1 


1 15.355 ±0.010 

2 15.230 ±0.050 

3 17.456 ± 0.043 

4 14.488 ± 0.018 

5 16.477 ±0.018 


14.738 ±0.002 
15.109 ±0.003 
16.788 ±0.007 
14.698 ± 0.002 
16.539 ±0.008 


13.885 ± 0.002 
14.447 ± 0.002 
15.918 ± 0.003 
14.334 ±0.002 
15.794 ± 0.003 


13.493 ± 0.002 
14.069 ± 0.002 
15.437 ± 0.005 
14.067 ± 0.002 
15.448 ± 0.005 


12.973 ± 0.003 
13.703 ±0.003 
15.013 ±0.007 
13.797 ±0.003 
15.125 ±0.008 
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Fig. 2.— Multi-band (UBVRI) light curves of 1ES 1011+496 from our MDM observations 
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Fig. 3. — Multi-band (UBVRI) light curves of 1ES 1011+496 during the high activity state 
in April - May 2008 
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However, we did observe substantial variability during our two observing runs in 2008 
April 9-14 and May 10 - 13. This was during the rising phase of a major outburst that 
peaked later that year. Unfortunately, to our knowledge, there are no observations available 
covering the peak of that outburst. During our observations we found a maximum R-band 
brightness of R ~ 14.8. The object later exceeded -R pea k < 14. q§. Fig. [3] shows the multi- 
band light curves from our 2008 April and May runs. They exhibit variability on time scales 
of a few days, but no evidence for intraday variability. The variability in all optical bands 
appears well correlated. This correlation will be investigated in more detail in §4J 



3. Optical spectral variability 

In order to test whether the variability discussed in the previous section is associated 
with spectral changes, we first calculated B - R color indices for any pair of B and R mag- 
nitudes measured within 15 minutes of each other. The resulting color-magnitude diagram 
(R magnitude vs. B - R color) is shown in Figure HI Error bars on B - R are calculated via 
standard error progatation, i.e., <Jb-r = ^/&r + The data clearly indicate color variabil- 
ity. A fit of a constant B - R color as a function of R magnitude results in xt = 2.86. A fit 
of a linear correlation results in a marginally acceptable xl = 1-36. A correlation analysis of 
the color-magnitude data set yields a Pearson's correlation coefficient of r = 0.57, which is 
generally interpreted as a weak positive correlation between color and magnitude. In order 
to quantify the probability of such a correlation coefficient resulting from an uncorrelated 
data set, we performed Monte-Carlo simulations of 1 billion randomly produced, uncorre- 
lated data sets, extending over similar spreads of values, and with the same number of data 
points as our observational data set. The Pearson's correlation coefficient for each set was 
evaluated, and from the entire ensemble, the probability of a correlation coefficient |r| > x 
for values of < x < 1 resulting from an uncorrelated data set was evaluated. From these 
simulations, we find that the probability of a correlation coefficient of r = 0.57 in a data set 
with the characteristics of our R vs. B - R data is P(\r\ > 0.57) ~ 4 x 10~ 9 , indicating, in 
fact, a highly significant correlation. The observed correlation corresponds to a bluer- when- 
brighter trend, as observed in most BL Lac objects. This is likely to reflect the dynamics of 
the non-thermal synchrotron emission from the jet dominating in the optical regime. 

In order to investigate spectral changes in the optical continuum in more detail, we 
extracted UBVRI SEDs for all sequences of magnitudes measured within 15 minutes of each 
other. For this purpose, the magnitudes were de-reddened using the Galactic extinction 



http : //users . utu. f i/kani/lm/ index .html 
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Fig. 4. — Color- magnitude diagram for 1ES 1011+496. The data show significant color 
variability. A positive linear correlation (bluer when brighter) is indicated by a Pearson's 
correlation coefficient of r = 0.57, with a probability for non-correlation of P(> r) ps 4x 10~ 9 . 
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Fig. 5.— Snap-shot UBVRI SEDs of 1ES 1011+496. The SEDs reveal a vF u peak typically 
between the B and V band and suggest a positive correlation between vF v peak flux and 
peak frequency. 
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coefficients as given in the NASA Extragalactic Databasej and converted to vF v fluxes. A 
representative sample of the resulting SEDs is plotted in Figure [5j The SEDs all exhibit 
a vF v peak in the optical regime, typically between the V and B bands. They suggest a 
positive trend of increasing vF v peak flux with increasing peak frequency, in accordance with 
the weak B - R vs. R correlation found above. We tested this hypothesis further by fitting 
all optical SEDs with a simple parabolic shape to determine the peak frequency, z/ pea k, and 
the peak flux, uF^ k . The best fit values for our entire data set are plotted in Fig. [6j 

The vF v peak flux and peak frequencies are clearly correlated, with Pearson's r = 0.78, 
with a probability for non-correlation of P(> r) w 10~ 15 . The best linear regression fit to 
the logarithms of the uF^ k and z/ p k values yields a power-law correlation vF^ k oc with 
k = 2.05 ± 0.17. A possible interpretation of this synchrotron peak shift will be discussed in 

m 

A visual inspection of Figure [6] seems to suggest a steepening of the peak flux vs. peak 
frequency dependence towards high peak frequencies (and peak fluxes). However, when 
restricting the regression to high frequencies (e.g., z/ p k > 6.2 x 10 14 Hz), the correlation 
between peak flux and peak frequency vanishes, therefore not allowing for a quantification 
of a possible change of the correlation slope towards high frequencies. 



4. Cross-correlation analysis 

A visual inspection of the light curves in Fig. [2] and [3] suggests that the variability in 
all wavebands is closely correlated. In order to corroborate this finding, we performed a 



Discrete Correlation Function (DCF, lEdelson fc Krolild Il988l ) analysis among all our light 



curves. Figure [7] shows a typical example of the resulting DCF between the B and R band 
light curves. The DCFs between the light curves of all bands peak at values near 1, indicating 
a close correlation between all optical bands. This is expected if the optical continuum is 
dominated by nonthermal synchrotron emission of the same relativistic electron population. 
Time lags between different frequency bands could potentially serve as a diagnos tic of the 



magn etic field strength in the emitting region (modulo the Doppler factor, see, e.g.. lBottcher 



20071 ). We therefore fitted the DCFs with an asymmetric Gaussian to determine possible 
inter-band delays through the fitted peak of the DCF. However, any lags indicated by our 
DCFs are all either consistent with or at the ~ 1 - 2 a level. Furthermore, our sequential 
data-taking process introduces an artificial "lag" of up to ~ 10 min, and none of the lags 
found through the DCF analysis are larger than that. Therefore, we conclude that we did 



http : //nedwww. ipac . caltech.edu/ 
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Fig. 6. — Best-fit peak frequency vs. peak vF v flux for the entire data set. The data show a 
significant correlation (r = 0.78, P(> r) ps 10~ 15 ), best fit by a dependence vF^ oc ^ k with 
k = 2.05 ± 0.17. The dashed line indicates a putative correlation with index k = 2. 
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Fig. 7. — Discrete Correlation Function between the B and R bands. A positive r would 
indicate a hard lag. The DCF has been fitted with an asymmetric Gaussian. The best-fit 
peak delay is r pk = (—4.8 ± 3.2) min. 



-14- 



not detect any inter-band time lags. 



5. Summary and Discussion 

We have presented an analysis of data from 5 years of observations of the BL Lac object 
1ES 1011+496 at the 1.3m McGraw-Hill Telescope of the MDM Observatory. We found 
moderate variability on a time scale of several days throughout most of our observations. 
The variability at all (UBVRI) optical bands is well correlated with no detectable time lags 
between them. The snap-shot SEDs during our observations showed a synchrotron peak 
within the optical range, typically between the V and B bands. The B - R color is correlated 
(r = 0.57; P U ncorr(> r) ~ 4 x 10 -8 ) with the R-band magnitude, indicating a bluer-when- 
brighter trend. Such a trend is observed in many BL Lac objects, where the optical emission 
is strongly dominated by synchrotron emission from the jet. We note that the opposite 
behaviour has been found in several quasar-type blazars, where a slowly variable Big Blue 
Bump, signaling a contribution due to a luminous acc retion disk, dilutes the continuum 



variability at the blue end of the optical spectrum (e.g., iRaiteri et all 120081 ). 



An analysis of the location of the synchrotron peak within the optical regime reveals a 
peak shift characterized by vF^ k oc z/p k with k = 2.05 ± 0.17, consistent with a z^ k scaling. 

There is a range of possible causes for the optical (and multi-wavelength) variability 
of blazar emission. These include changes in the Doppler factor (e.g., caused by a bending 
jet), injection of a new relativistic particle population into the jet (plausibly caused by 
a shock), a changing acceleration efficiency (changing characteristic Lorentz factors of the 
radiating electrons and possibly a change of the spectral index of the non-thermal electron 
distribution), and/or a change of the magnetic field. In a realistic scenario, several of these 
effects might be at work at the same time to produce the observed blazar variability. However, 
one can make simple predictions concerning the shift (in frequency and vF v peak flux) of the 
synchrotron peak for at least three cases: A changing Doppler factor, a changing magnetic 
field, and a change of the characteristic (peak) electron Lorentz factor (leaving all other 
parameters of the emission region unchanged). 

The peak frequency of the synchrotron spectrum is related to a peak in the electron 
spectrum at a characteristic Lorentz factor j p , the magnetic field B and the Doppler factor 
D = (r [1 - f3 T cos flobs])" 1 , where T = (1 - /^)~ 1/2 is the bulk Lorentz factor of the emission 
region and ^obs is the angle between the line of sight towards Earth and the direction of 
motion (the jet axis), through 
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p pk cxY p BD (1) 
The vF v flux at the synchrotron peak is related to those quantities through 

uFf oc 7p 2 B 2 D 4 (2) 

From equations [1] and |2], we see that if the variability is dominated by a changing Doppler 
factor, one would expect a synchrotron peak shift as vF^ oc v 4 k . A change solely in the 
characteristic electron Lorentz factor, j p , would result in a synchrotron peak shift as uF^ k oc 
Upk, while a change in only the magnetic field yields the behaviour uF^ k oc z^ k . 

Therefore, we conclude that the synchrotron peak shift found in our data set is consistent 
with the variability being dominated by a changing magnetic field. However, as pointed out 
above, we need to caution that such a change in the magnetic field might realistically also 
impact the shape of the electron distribution, primarily through a changing synchrotron 
cooling time scale. Clearly, more sophisticated analyses of this synchrotron peak shift are 
needed, but are beyond the scope of this paper. 

We note that the shape of the optical SEDs found in all of our observations contradicts 
the optical-UV spectrum observed by Swift/U VOT on 2008 May 2 and 8, which indicates 



a rising slope throughout the optical regime (lAbdo et al.l 120101 ). However, most of our 
observations were taken during moderately faint states of the source, while the Swift/UVOT 
spectrum corresponds to a bright state, similar to the major optical flare that triggered the 
MAGIC detection in 2007. Given the trend of the synchrotron peak shift which we found 
in our data, it is conceivable that the Swift/UVOT observations correspond to an extreme 
case of a high synchrotron peak frequency, in accord with a very high optical flux. 
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